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1 Project summary 

This research aimed at applying calculating the sound emission from a driven cavity using 
LES. 

2 Accomplished tasks 

During the course of this research, first we have determined the accuracy of the available 
subgrid-scale models in non-equilibrium flows. It was found that the Lagrangian eddy- 
viscosity model and the Lagrangian mixed models give the most accurate results. These 
results were published in the two enclosed articles [2, 1], 

Secondly, we attempted to use a spectral multi-domain code for the problem under con- 
sideration. The effort was, however, unsuccessful. The code developed instabilities at the 
sharp corners of the cavity due to discontinuities in the velocity derivative. 

Then, a finite-difference, staggered code was developed. The code was tested by comput- 
ing the flow-field in the flat-plate boundary layer studied numerically and experimentally by 
Spalart and Watmuff [3]. The immersed-boundary technique [4, 5] has been implemented 
and tested both in laminar and turbulent cases. This technique consists in applying dis- 
tributed body forces inside the computational domain. These forces are designed in such a 
way that the fluid velocity goes to zero (or to some assigned value) at specified points. 

The “Direct forcing” method [5] has been used. This method consists in the following 
steps: within the fractional time-step framework, the velocity is solved from the Helmholtz 
equation: 
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The body force is then assigned to be zero everywhere except on the immersed-body surface, 
T, where it is given by 

fiM = - RHS. + ^p, 

where V t is the body velocity. This forces the no-slip condition to hold on the predicted 
velocity at the internal points. The body force was implemented va a linear interpolation of 
the velocity between the body itself and the second grid-point outside the body[5]. 

With this method, the pressure gradient on the immersed body itself is zero; this can be 
shown very easily by considering the Navier-Stokes equations in a r - n frame of reference 
that moves with the body, where (r is tangent, n normal to the body surface): 
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It is easy to show that all the terms on the left-hand-side of the equation, as well as 
the tangential viscous term are identically zero. If the velocity is obtained from linear 
interpolation in the direction normal to the body, as is the case in the present implementation, 
the remaining viscous term is also zero, and the pressure satisfies Vp = 0 on T. 

A configuration close to the problem of interest, a backward-facing step, has been cal- 
culated. The results appear to be in good agreement with the experimental data. The 
reattachment point, for instance (Fig. 1) is predicted correctly. 

Unfortunately, due to the significant amount of time spent dealing with the problems with 
the spectral code, the proposed cavity calculation could not be completed. The development 
of the present code, however, will be beneficial for future work of interest to NASA. 
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